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In light of the mounting evidence that the highest energy cosmic rays are dominated by protons and not gamma- 
rays, we discuss the prospect that these cosmic rays are generated in the decay or annihilation of superheavy relic 
particles. We calculate the high energy neutrino spectrum which results and normalize our results to the ultra- 
high energy cosmic ray spectrum. We show that most sceneries are already constrained by present limits placed 
by the AMANDA experiment. 



1. Introduction 

The discovery of cosmic rays with energy ex- 
ceeding the GZK cutoff presents an interest- 
ing challenge to astrophysics, particle physics, or 
both. Numerous scenarios have been proposed to 
solve the problem. These include exotic particles 
j2|, neutrinos with QCD scale cross sections 
semi-local astrophysical sources [Q and top-down 
models 

Recent measurements confirm that our uni- 
verse contains a large fraction of cold dark mat- 
ter 1^]. A top-down model in which annihilat- 
ing or decaying superheavy particles produce the 
highest energy cosmic rays could potentially solve 
both of these problems . 

Conventional particle physics implies that ul- 
tra high-energy jets fragment predominantly into 
photons with a small admixture of protons. This 
seems to be in disagreement with mounting evi- 
dence that the highest energy cosmic rays are not 
photons 1^. This does not necessarily rule out 
superheavy particles as the source of the highest 
energy cosmic rays. The uncertainties associated 
with the cascading of the jets in the universal 
radio background and with the strength of inter- 
galactic magnetic fields leave open the possibility 
that ultra high-energy photons may be depleted 
from the cosmic ray spectrum near 10^*^ eV, leav- 
ing a dominant proton component at GZK en- 
ergies [^. With this in mind, we will choose 
to normalize the proton spectrum from top-down 
scenarios with the observed ultra high-energy cos- 



mic ray flux. 

Neutrinos are produced more numerously than 
protons and travel much greater distances. The 
main point of this note is to point out that this 
"renormalization" of the observed cosmic ray flux 
to protons generically predicts observable neu- 
trino signals in operating high-energy neutrino 
telescopes such as AMANDA. 

2. Nucleons from Ultra High-Energy Jets 

To normalize the production rate of ultra high- 
energy jets, it is necessary to calculate the spec- 
trum of nucleons resulting from their fragmenta- 
tion. Each jet will fragment into a large num- 
ber of hadrons approximated by a fragmenta- 
tion function rooted in accelerator data ||ll|. All 
hadrons produced eventually decay into pions and 
nucleons. For a detailed discussion of ultra high- 
energy fragmentation, see Ref. p2[ . 

To solve the ultra high-energy cosmic ray 
problem, this nucleon flux must accommodate 
the events above the GZK cutoff. Observa- 
tions indicate on the order of 10^^^ events 
m~^s~^sr" -iGeV"^ in the energy range above the 
GZK cutoff (5 X 1019 eV to 2 X 10^° eV)|l|. The 
formalism of a generic top-down scenario is suffi- 
ciently flexible to explain the data from either the 
HIRES or AGASA experiments. For an example, 
see Figure 1. 
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Figure 1. The ultra high-energy cosmic ray flux 
predicted from the decay or annihilation of super- 
heavy dark matter particles producing 10^^ eV 
hadronic jets is compared to the AGASA cosmic 
ray data. The distribution of dark matter used is 
isotropic with an overdensity factor of 10^ within 
20 kpc. Note that all observed super GZK events 
can be explained by this mechanism. 

3. Cross Sections and Lifetimes for Super- 
heavy Particles as Dark Matter 

If the superheavy relics responsible for the 
highest energy cosmic rays are also the solution 
to the dark matter problem, then 

PX ^ 0.3 X peritical ^ 0.3 X 8.45 X lO'^kg/m^, ( 

and may be much larger locally. Therefore, the 
lifetime of such a decaying particle must be: 

PX 



TX 



10^'' years. 



1) 



(2) 



Such a long lifetime may be disfavored by fine- 
tuning arguments, but can be possible If, 
however, the superheavy relics in question were 
not the major dark matter component, this life- 
time could be much shorter. 

If instead we consider stable particles which can 
annihilate with each other, we can calculate their 
annihilation cross section as a function of velocity. 
For an isotropic distribution of particles, 

dux 



dt 



f^xx 



2 X lO-i^m^s^i 



or, for a galactic halo distribution. 



c^xx 



lO-i^m^s-i 



(4) 



(5) 



A characteristic velocity of 500 km/s, for exam- 
ple, would correspond to a annihilation cross sec- 
tion of ^ 10^ bn for extragalactic dark matter or 
^ 10"^ bn for a galactic halo distribution. If the 
dark matter were not uniformly distributed, how- 
ever, but were distributed in clumps with char- 
acteristic densities of /Ociump Cpmcan then the 
annihilation cross section would be lowered by 
a factor of C. For example, if the density of a 
typical cluster were 10® times greater than in an 
isotropic distribution, annihilation cross sections 
could be in the mb range. It is interesting to 
note that if superheavy dark matter is distributed 
locally with mb elastic scattering cross sections, 
they would become gravitationally trapped in the 
sun and annihilate as described in Rcf. [isl] . 

4. Neutrinos from Ultra High-Energy Jets 

Neutrinos are produced in several ways in 
the fragmentation of ultra high-energy jets: in 
semileptonic bottom and charm decays, in W pro- 
ducing top decays and most importantly, in the 
decay of charged pions. For a detailed discussion, 
see Ref. (l|]. 

To obtain the neutrino flux, we multiply the 
injection spectrum by the average distance trav- 
eled by a neutrino and by the rate per volume for 
hadronic jets which we calculated earlier. Neu- 
trinos, not being limited by scattering, travel up 
to the age of the universe at the speed of light 
(~ 3000 Mpc in an Euclidean approximation). 
The predicted neutrino flux is shown in Figure 
2. Note that this flux is significantly higher than 
the present limits placed by the AMANDA exper- 
iment for the case of an isotropic distribution 
of ultra high-energy jets. For the galactic sce- 
nario, however, the flux predicted is comparable 
with present AMANDA-BIO hmits. Also shown 
are the limits anticipated from AMANDA-II data 
and the IceCube experiment. 
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Figure 2. The muon neutrino flux from 10^^ eV 
jets normalized to the highest energy cosmic rays. 
Shown are the fluxes calculated for an isotropic 
and for a galactic distribution of jets. The dotted 
lines shown are the experimental diffuse flux lim- 
its from present AMANDA-BIO data, projected 
AMANDA-II data and projected IceCube data. 
These limits are < 9 x 10"'^, 9 x 10"^ and 

5 X 10~^ GeV cm~^s~^sr~^ respectively. Note 
that operating experiments can effectively test 
this range of models. 



5. Event Rates in High-energy Neutrino 
Telescopes 

The diffuse flux of high-energy neutrinos can 
be observed by operating and planned neutrino 
telescopes. AMANDA-BIO, with an effective area 
of ~5,000 square meters has placed the strongest 
limits on the diffuse neutrino flux. Figure 2 shows 
that for a galactic distribution of superheavy dark 
matter particles, the predicted high-energy neu- 
trino flux is on the order of the present diffuse 
flux limit. The neutrino flux for an isotropic dis- 
tribution of ultra high-energy jets is well above 
the present limits. If superheavy particles are 
indeed the source of the highest energy cosmic 
rays, observations by high-energy neutrino tele- 
scopes should reveal the corresponding neutrino 
flux. For a review of cerenkov neutrino telescopes 



seeRef. 

Note that gamma ray astronomy can also pro- 
vide an interesting test of these scenarios . 

6. Conclusions 

If the decay or annihilation of superheavy relics 
is the source of the highest energy cosmic rays, 
then a high-energy neutrino flux should accom- 
pany the observed cosmic ray flux. This neutrino 
flux will be much higher than the flux of nucleons 
due to the much greater mean free path of neu- 
trinos and greater multiplicity of neutrinos pro- 
duced in high-energy hadronic jets. 

The high-energy neutrino flux generated in 
such a scenario can be calculated by normaliz- 
ing the flux of appropriate particles to the ultra 
high-energy cosmic ray flux. With mounting evi- 
dence that the highest energy cosmic rays are pro- 
tons or nuclei and not photons, we must require 
that the ultra high-energy photons are degraded 
by the universal radio background, leaving pro- 
tons to dominate the highest energy cosmic ray 
flux. The neutrino flux must be normalized to 
the proton flux resulting in significantly improved 
prospects for its detection. 

This paper shows that the neutrino flux accom- 
panying the highest energy cosmic rays in these 
models is on the order of the limits placed by op- 
erating high-energy neutrino telescopes such as 
AMANDA. Further data from AMANDA II, or 
next generation neutrino telescope IceCube, can 
test the viability of models in which superheavy 
particles generate the highest energy cosmic rays. 
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